Genomics Research of Lifetime Depression in the Netherlands: The BIObanks Netherlands Internet Collaboration (BIONIC) Project.

In this cohort profile article we describe the lifetime major depressive disorder (MDD) database that has been established as part of the BIObanks Netherlands Internet Collaboration (BIONIC). Across the Netherlands we collected data on Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition (DSM-5) lifetime MDD diagnosis in 132,850 Dutch individuals. Currently, N = 66,684 of these also have genomewide single nucleotide polymorphism (SNP) data. We initiated this project because the complex genetic basis of MDD requires large population-wide studies with uniform in-depth phenotyping. For standardized phenotyping we developed the LIDAS (LIfetime Depression Assessment Survey), which then was used to measure MDD in 11 Dutch cohorts. Data from these cohorts were combined with diagnostic interview depression data from 5 clinical cohorts to create a dataset of N = 29,650 lifetime MDD cases (22%) meeting DSM-5 criteria and 94,300 screened controls. In addition, genomewide genotype data from the cohorts were assembled into a genomewide association study (GWAS) dataset of N = 66,684 Dutch individuals (25.3% cases). Phenotype data include DSM-5-based MDD diagnoses, sociodemographic variables, information on lifestyle and BMI, characteristics of depressive symptoms and episodes, and psychiatric diagnosis and treatment history. We describe the establishment and harmonization of the BIONIC phenotype and GWAS datasets and provide an overview of the available information and sample characteristics. Our next step is the GWAS of lifetime MDD in the Netherlands, with future plans including fine-grained genetic analyses of depression characteristics, international collaborations and multi-omics studies.

Association between circulating inflammatory markers and adult cancer risk: a Mendelian randomization analysis.

BACKGROUND: Tumour-promoting inflammation is a "hallmark" of cancer and conventional epidemiological studies have reported links between various inflammatory markers and cancer risk. The causal nature of these relationships and, thus, the suitability of these markers as intervention targets for cancer prevention is unclear. METHODS: We meta-analysed 6 genome-wide association studies of circulating inflammatory markers comprising 59,969 participants of European ancestry. We then used combined cis-Mendelian randomization and colocalisation analysis to evaluate the causal role of 66 circulating inflammatory markers in risk of 30 adult cancers in 338,294 cancer cases and up to 1,238,345 controls. Genetic instruments for inflammatory markers were constructed using genome-wide significant (P < 5.0 × 10-8) cis-acting SNPs (i.e., in or ±250 kb from the gene encoding the relevant protein) in weak linkage disequilibrium (LD, r2 < 0.10). Effect estimates were generated using inverse-variance weighted random-effects models and standard errors were inflated to account for weak LD between variants with reference to the 1000 Genomes Phase 3 CEU panel. A false discovery rate (FDR)-corrected P-value ("q-value") <0.05 was used as a threshold to define "strong evidence" to support associations and 0.05 ≤ q-value < 0.20 to define "suggestive evidence". A colocalisation posterior probability (PPH4) >70% was employed to indicate support for shared causal variants across inflammatory markers and cancer outcomes. Findings were replicated in the FinnGen study and then pooled using meta-analysis. FINDINGS: We found strong evidence to support an association of genetically-proxied circulating pro-adrenomedullin concentrations with increased breast cancer risk (OR: 1.19, 95% CI: 1.10-1.29, q-value = 0.033, PPH4 = 84.3%) and suggestive evidence to support associations of interleukin-23 receptor concentrations with increased pancreatic cancer risk (OR: 1.42, 95% CI: 1.20-1.69, q-value = 0.055, PPH4 = 73.9%), prothrombin concentrations with decreased basal cell carcinoma risk (OR: 0.66, 95% CI: 0.53-0.81, q-value = 0.067, PPH4 = 81.8%), and interleukin-1 receptor-like 1 concentrations with decreased triple-negative breast cancer risk (OR: 0.92, 95% CI: 0.88-0.97, q-value = 0.15, PPH4 = 85.6%). These findings were replicated in pooled analyses with the FinnGen study. Though suggestive evidence was found to support an association of macrophage migration inhibitory factor concentrations with increased bladder cancer risk (OR: 2.46, 95% CI: 1.48-4.10, q-value = 0.072, PPH4 = 76.1%), this finding was not replicated when pooled with the FinnGen study. For 22 of 30 cancer outcomes examined, there was little evidence (q-value ≥0.20) that any of the 66 circulating inflammatory markers examined were associated with cancer risk. INTERPRETATION: Our comprehensive joint Mendelian randomization and colocalisation analysis of the role of circulating inflammatory markers in cancer risk identified potential roles for 4 circulating inflammatory markers in risk of 4 site-specific cancers. Contrary to reports from some prior conventional epidemiological studies, we found little evidence of association of circulating inflammatory markers with the majority of site-specific cancers evaluated. FUNDING: Cancer Research UK (C68933/A28534, C18281/A29019, PPRCPJT∖100005), World Cancer Research Fund (IIG_FULL_2020_022), National Institute for Health Research (NIHR202411, BRC-1215-20011), Medical Research Council (MC_UU_00011/1, MC_UU_00011/3, MC_UU_00011/6, and MC_UU_00011/4), Academy of Finland Project 326291, European Union's Horizon 2020 grant agreement no. 848158 (EarlyCause), French National Cancer Institute (INCa SHSESP20, 2020-076), Versus Arthritis (21173, 21754, 21755), National Institutes of Health (U19 CA203654), National Cancer Institute (U19CA203654).

Multi-trait analysis characterizes the genetics of thyroid function and identifies causal associations with clinical implications.

To date only a fraction of the genetic footprint of thyroid function has been clarified. We report a genome-wide association study meta-analysis of thyroid function in up to 271,040 individuals of European ancestry, including reference range thyrotropin (TSH), free thyroxine (FT4), free and total triiodothyronine (T3), proxies for metabolism (T3/FT4 ratio) as well as dichotomized high and low TSH levels. We revealed 259 independent significant associations for TSH (61% novel), 85 for FT4 (67% novel), and 62 novel signals for the T3 related traits. The loci explained 14.1%, 6.0%, 9.5% and 1.1% of the total variation in TSH, FT4, total T3 and free T3 concentrations, respectively. Genetic correlations indicate that TSH associated loci reflect the thyroid function determined by free T3, whereas the FT4 associations represent the thyroid hormone metabolism. Polygenic risk score and Mendelian randomization analyses showed the effects of genetically determined variation in thyroid function on various clinical outcomes, including cardiovascular risk factors and diseases, autoimmune diseases, and cancer. In conclusion, our results improve the understanding of thyroid hormone physiology and highlight the pleiotropic effects of thyroid function on various diseases.

Causal relationships between risk of venous thromboembolism and 18 cancers: a bidirectional Mendelian randomization analysis.

BACKGROUND: People with cancer experience high rates of venous thromboembolism (VTE). Risk of subsequent cancer is also increased in people experiencing their first VTE. The causal mechanisms underlying this association are not completely understood, and it is unknown whether VTE is itself a risk factor for cancer. METHODS: We used data from large genome-wide association study meta-analyses to perform bidirectional Mendelian randomization analyses to estimate causal associations between genetic liability to VTE and risk of 18 different cancers. RESULTS: We found no conclusive evidence that genetic liability to VTE was causally associated with an increased incidence of cancer, or vice versa. We observed an association between liability to VTE and pancreatic cancer risk [odds ratio for pancreatic cancer: 1.23 (95% confidence interval: 1.08-1.40) per log-odds increase in VTE risk, P = 0.002]. However, sensitivity analyses revealed this association was predominantly driven by a variant proxying non-O blood group, with inadequate evidence to suggest a causal relationship. CONCLUSIONS: These findings do not support the hypothesis that genetic liability to VTE is a cause of cancer. Existing observational epidemiological associations between VTE and cancer are therefore more likely to be driven by pathophysiological changes which occur in the setting of active cancer and anti-cancer treatments. Further work is required to explore and synthesize evidence for these mechanisms.

Whole Genome Sequencing Identifies Novel Common and Low-Frequency Variants Associated With Age-Related Macular Degeneration.

Purpose: To identify associations of common, low-frequency, and rare variants with advanced age-related macular degeneration (AMD) using whole genome sequencing (WGS). Methods: WGS data were obtained for 2123 advanced AMD patients (participants of clinical trials for advanced AMD) and 2704 controls (participants of clinical trials for asthma [N = 2518] and Alzheimer's disease [N = 186]), and joint genotype calling was performed, followed by quality control of the dataset. Single variant association analyses were performed for all identified common, low-frequency, and rare variants. Gene-based tests were executed for rare and low-frequency variants using SKAT-O and three groups of variants based on putative impact information: (1) all variants, (2) modifier impact variants, and (3) high- and moderate-impact variants. To ascertain independence of the identified associations from previously reported AMD and asthma loci, conditional analyses were performed. Results: Previously identified AMD variants at the CFH, ARMS2/HTRA1, APOE, and C3 loci were associated with AMD at a genome-wide significance level. We identified new single variant associations for common variants near the PARK7 gene and in the long non-coding RNA AC103876.1, and for a rare variant near the TENM3 gene. In addition, gene-based association analyses identified a burden of modifier variants in eight intergenic and gene-spanning regions and of high- and moderate-impact variants in the C3, CFHR5, SLC16A8, and CFI genes. Conclusions: We describe the largest WGS study in AMD to date. We confirmed previously identified associations and identified several novel associations that are worth exploring in further follow-up studies.

Genetic insights into resting heart rate and its role in cardiovascular disease.

Resting heart rate is associated with cardiovascular diseases and mortality in observational and Mendelian randomization studies. The aims of this study are to extend the number of resting heart rate associated genetic variants and to obtain further insights in resting heart rate biology and its clinical consequences. A genome-wide meta-analysis of 100 studies in up to 835,465 individuals reveals 493 independent genetic variants in 352 loci, including 68 genetic variants outside previously identified resting heart rate associated loci. We prioritize 670 genes and in silico annotations point to their enrichment in cardiomyocytes and provide insights in their ECG signature. Two-sample Mendelian randomization analyses indicate that higher genetically predicted resting heart rate increases risk of dilated cardiomyopathy, but decreases risk of developing atrial fibrillation, ischemic stroke, and cardio-embolic stroke. We do not find evidence for a linear or non-linear genetic association between resting heart rate and all-cause mortality in contrast to our previous Mendelian randomization study. Systematic alteration of key differences between the current and previous Mendelian randomization study indicates that the most likely cause of the discrepancy between these studies arises from false positive findings in previous one-sample MR analyses caused by weak-instrument bias at lower P-value thresholds. The results extend our understanding of resting heart rate biology and give additional insights in its role in cardiovascular disease development.

Causal relationships between risk of venous thromboembolism and 18 cancers: a bidirectional Mendelian randomisation analysis.

Background: People with cancer experience high rates of venous thromboembolism (VTE). Additionally, risk of subsequent cancer is increased in people experiencing their first VTE. The causal mechanisms underlying this association are not completely understood, and it is unknown whether VTE is itself a risk factor for cancer. Methods: We used data from large genome-wide association study meta-analyses to perform bi-directional Mendelian randomisation analyses to estimate causal associations between genetically-proxied lifetime risk of VTE and risk of 18 different cancers. Results: We found no conclusive evidence that genetically-proxied lifetime risk of VTE was causally associated with an increased incidence of cancer, or vice-versa. We observed an association between VTE and pancreatic cancer risk (odds ratio for pancreatic cancer 1.23 (95% confidence interval 1.08 - 1.40) per log-odds increase in risk of VTE, P = 0.002). However, sensitivity analyses revealed this association was predominantly driven by a variant proxying non-O blood group, with inadequate evidence from Mendelian randomisation to suggest a causal relationship. Conclusions: These findings do not support the hypothesis that genetically-proxied lifetime risk of VTE is a cause of cancer. Existing observational epidemiological associations between VTE and cancer are therefore more likely to be driven by pathophysiological changes which occur in the setting of active cancer and anti-cancer treatments. Further work is required to explore and synthesise evidence for these mechanisms.

Association between circulating inflammatory markers and adult cancer risk: a Mendelian randomization analysis.

Background: Tumour-promoting inflammation is a "hallmark" of cancer and conventional epidemiological studies have reported links between various inflammatory markers and cancer risk. The causal nature of these relationships and, thus, the suitability of these markers as intervention targets for cancer prevention is unclear. Methods: We meta-analysed 6 genome-wide association studies of circulating inflammatory markers comprising 59,969 participants of European ancestry. We then used combined cis-Mendelian randomization and colocalisation analysis to evaluate the causal role of 66 circulating inflammatory markers in risk of 30 adult cancers in 338,162 cancer cases and up to 824,556 controls. Genetic instruments for inflammatory markers were constructed using genome-wide significant (P < 5.0 x 10-8) cis-acting SNPs (i.e. in or ±250 kb from the gene encoding the relevant protein) in weak linkage disequilibrium (LD, r2 < 0.10). Effect estimates were generated using inverse-variance weighted random-effects models and standard errors were inflated to account for weak LD between variants with reference to the 1000 Genomes Phase 3 CEU panel. A false discovery rate (FDR)-corrected P-value ("q-value") < 0.05 was used as a threshold to define "strong evidence" to support associations and 0.05 ≤ q-value < 0.20 to define "suggestive evidence". A colocalisation posterior probability (PPH4) > 70% was employed to indicate support for shared causal variants across inflammatory markers and cancer outcomes. Results: We found strong evidence to support an association of genetically-proxied circulating pro-adrenomedullin concentrations with increased breast cancer risk (OR 1.19, 95% CI 1.10-1.29, q-value=0.033, PPH4=84.3%) and suggestive evidence to support associations of interleukin-23 receptor concentrations with increased pancreatic cancer risk (OR 1.42, 95% CI 1.20-1.69, q-value=0.055, PPH4=73.9%), prothrombin concentrations with decreased basal cell carcinoma risk (OR 0.66, 95% CI 0.53-0.81, q-value=0.067, PPH4=81.8%), macrophage migration inhibitory factor concentrations with increased bladder cancer risk (OR 1.14, 95% CI 1.05-1.23, q-value=0.072, PPH4=76.1%), and interleukin-1 receptor-like 1 concentrations with decreased triple-negative breast cancer risk (OR 0.92, 95% CI 0.88-0.97, q-value=0.15), PPH4=85.6%). For 22 of 30 cancer outcomes examined, there was little evidence (q-value ≥ 0.20) that any of the 66 circulating inflammatory markers examined were associated with cancer risk. Conclusion: Our comprehensive joint Mendelian randomization and colocalisation analysis of the role of circulating inflammatory markers in cancer risk identified potential roles for 5 circulating inflammatory markers in risk of 5 site-specific cancers. Contrary to reports from some prior conventional epidemiological studies, we found little evidence of association of circulating inflammatory markers with the majority of site-specific cancers evaluated.

Fast and accurate recurrent event analysis for genome-wide association studies.

Many diseases recur after recovery, for example, recurrences in cancer and infections. However, research is often focused on analysing only time-to-first recurrence, thereby ignoring any subsequent recurrences that may occur after the first. Statistical models for the analysis of recurrent events are available, of which the extended Cox proportional hazards frailty model is the current state-of-the-art. However, this model is too statistically complex for computationally efficient application in high-dimensional data sets, including genome-wide association studies (GWAS). Here, we develop an application for fast and accurate recurrent event analysis in GWAS, called SPARE (SaddlePoint Approximation for Recurrent Event analysis). In SPARE, every DNA variant is tested for association with recurrence risk using a modified score statistic. A saddlepoint approximation is implemented to achieve statistical accuracy. SPARE controls the Type I error, and its statistical power is similar to existing recurrent event models, yet SPARE is significantly faster. An application of SPARE in a recurrent event GWAS on bladder cancer for 6.2 million DNA variants in 1,443 individuals required less than 15 min, whereas existing recurrent event methods would require several weeks.

Intravesical BCG in patients with non-muscle invasive bladder cancer induces trained immunity and decreases respiratory infections.

BACKGROUND: BCG is recommended as intravesical immunotherapy to reduce the risk of tumor recurrence in patients with non-muscle invasive bladder cancer (NMIBC). Currently, it is unknown whether intravesical BCG application induces trained immunity. METHODS: The aim of this research was to determine whether BCG immunotherapy induces trained immunity in NMIBC patients. We conducted a prospective observational cohort study in 17 NMIBC patients scheduled for BCG therapy and measured trained immunity parameters at 9 time points before and during a 1-year BCG maintenance regimen. Ex vivo cytokine production by peripheral blood mononuclear cells, epigenetic modifications, and changes in the monocyte transcriptome were measured. The frequency of respiratory infections was investigated in two larger cohorts of BCG-treated and non-BCG treated NMIBC patients as a surrogate measurement of trained immunity. Gene-based association analysis of genetic variants in candidate trained immunity genes and their association with recurrence-free survival and progression-free survival after BCG therapy was performed to investigate the hypothesized link between trained immunity and clinical response. RESULTS: We found that intravesical BCG does induce trained immunity based on an increased production of TNF and IL-1ß after heterologous ex vivo stimulation of circulating monocytes 6-12 weeks after intravesical BCG treatment; and a 37% decreased risk (OR 0.63 (95% CI 0.40 to 1.01)) for respiratory infections in BCG-treated versus non-BCG-treated NMIBC patients. An epigenomics approach combining chromatin immuno precipitation-sequencing and RNA-sequencing with in vitro trained immunity experiments identified enhanced inflammasome activity in BCG-treated individuals. Finally, germline variation in genes that affect trained immunity was associated with recurrence and progression after BCG therapy in NMIBC. CONCLUSION: We conclude that BCG immunotherapy induces trained immunity in NMIBC patients and this may account for the protective effects against respiratory infections. The data of our gene-based association analysis suggest that a link between trained immunity and oncological outcome may exist. Future studies should further investigate how trained immunity affects the antitumor immune responses in BCG-treated NMIBC patients.

Sociodemographic, lifestyle and clinical characteristics of energy-related depression symptoms: A pooled analysis of 13,965 depressed cases in 8 Dutch cohorts.

BACKGROUND: In a substantial subgroup of depressed patients, atypical, energy-related depression symptoms (e.g. increased appetite/weight, hypersomnia, loss of energy) tend to cluster with immuno-metabolic dysregulations (e.g. increased BMI and inflammatory markers). This clustering is proposed to reflect a more homogeneous depression pathology. This study examines to what extent energy-related symptoms are associated and share sociodemographic, lifestyle and clinical characteristics. METHODS: Data were available from 13,965 participants from eight Dutch cohorts with DSM-5 lifetime major depression assessed by the Lifetime Depression Assessment Self-report (LIDAS) questionnaire. Information on four energy-related depression symptoms were extracted: energy loss, increased appetite, increased weight, and hypersomnia. Tetrachoric correlations between these symptoms, and associations of these symptoms with sociodemographic (sex, age, education), lifestyle (physical activity, BMI, smoking) and clinical characteristics (age of onset, episode duration, history, treatment and recency, and self-reported comorbidity) were computed. RESULTS: Correlations between energy-related symptoms were overall higher than those with other depression symptoms and varied from 0.90 (increased appetite vs increased weight) to 0.11 (increased appetite vs energy loss). All energy-related symptoms were strongly associated with higher BMI and a more severe clinical profile. Patients with increased appetite were more often smokers, and only patients with increased appetite or weight more often had a self-reported diagnosis of PTSD (OR = 1.17, p = 2.91E-08) and eating disorder (OR = 1.40, p = 4.08E-17). CONCLUSIONS: The symptom-specific associations may have consequences for a profile integrating these symptoms, which can be used to reflect immuno-metabolic depression. They indicate the need to study immuno-metabolic depression at individual symptom resolution as a starting point.

A Genome-Wide Association Study into the Aetiology of Congenital Solitary Functioning Kidney.

Congenital solitary functioning kidney (CSFK) is a birth defect that occurs in 1:1500 children and predisposes them to kidney injury. Its aetiology is likely multifactorial. In addition to known monogenic causes and environmental risk factors, common genetic variation may contribute to susceptibility to CSFK. We performed a genome-wide association study among 452 patients with CSFK and two control groups of 669 healthy children and 5363 unaffected adults. Variants in two loci reached the genome-wide significance threshold of 5 × 10-8, and variants in 30 loci reached the suggestive significance threshold of 1 × 10-5. Of these, an identified locus with lead single nucleotide variant (SNV) rs140804918 (odds ratio 3.1, p-value = 1.4 × 10-8) on chromosome 7 was most promising due to its close proximity to HGF, a gene known to be involved in kidney development. Based on their known molecular functions, both KCTD20 and STK38 could explain the suggestive significant association with lead SNV rs148413365 on chromosome 6. Our findings need replication in an independent cohort of CSFK patients before they can be established definitively. However, our analysis suggests that common variants play a role in CSFK aetiology. Future research could enhance our understanding of the molecular mechanisms involved.

The Effects of Common Genetic Variation in 96 Genes Involved in Thyroid Hormone Regulation on TSH and FT4 Concentrations.

OBJECTIVE: While most of the variation in thyroid function is determined by genetic factors, single nucleotide polymorphisms (SNPs) identified via genome-wide association analyses have only explained ~5% to 9% of this variance so far. Most SNPs were in or nearby genes with no known role in thyroid hormone (TH) regulation. Therefore, we performed a large-scale candidate gene study investigating the effect of common genetic variation in established TH regulating genes on serum thyrotropin [thyroid-stimulating hormone (TSH)] and thyroxine (FT4) concentrations. METHODS: SNPs in or within 10 kb of 96 TH regulating genes were included (30 031 TSH SNPs, and 29 962 FT4 SNPs). Associations were studied in 54 288 individuals from the ThyroidOmics Consortium. Linkage disequilibrium-based clumping was used to identify independently associated SNPs. SNP-based explained variances were calculated using SumHer software. RESULTS: We identified 23 novel TSH-associated SNPs in predominantly hypothalamic-pituitary-thyroid axis genes and 25 novel FT4-associated SNPs in mainly peripheral metabolism and transport genes. Genome-wide SNP variation explained ~21% (SD 1.7) of the total variation in both TSH and FT4 concentrations, whereas SNPs in the 96 TH regulating genes explained 1.9% to 2.6% (SD 0.4). CONCLUSION: Here we report the largest candidate gene analysis on thyroid function, resulting in a substantial increase in the number of genetic variants determining TSH and FT4 concentrations. Interestingly, these candidate gene SNPs explain only a minor part of the variation in TSH and FT4 concentrations, which substantiates the need for large genetic studies including common and rare variants to unravel novel, yet unknown, pathways in TH regulation.

Transferrin Saturation/Hepcidin Ratio Discriminates TMPRSS6-Related Iron Refractory Iron Deficiency Anemia from Patients with Multi-Causal Iron Deficiency Anemia.

Pathogenic TMPRSS6 variants impairing matriptase-2 function result in inappropriately high hepcidin levels relative to body iron status, leading to iron refractory iron deficiency anemia (IRIDA). As diagnosing IRIDA can be challenging due to its genotypical and phenotypical heterogeneity, we assessed the transferrin saturation (TSAT)/hepcidin ratio to distinguish IRIDA from multi-causal iron deficiency anemia (IDA). We included 20 IRIDA patients from a registry for rare inherited iron disorders and then enrolled 39 controls with IDA due to other causes. Plasma hepcidin-25 levels were measured by standardized isotope dilution mass spectrometry. IDA controls had not received iron therapy in the last 3 months and C-reactive protein levels were <10.0 mg/L. IRIDA patients had significantly lower TSAT/hepcidin ratios compared to IDA controls, median 0.6%/nM (interquartile range, IQR, 0.4-1.1%/nM) and 16.7%/nM (IQR, 12.0-24.0%/nM), respectively. The area under the curve for the TSAT/hepcidin ratio was 1.000 with 100% sensitivity and specificity (95% confidence intervals 84-100% and 91-100%, respectively) at an optimal cut-off point of 5.6%/nM. The TSAT/hepcidin ratio shows excellent performance in discriminating IRIDA from TMPRSS6-unrelated IDA early in the diagnostic work-up of IDA provided that recent iron therapy and moderate-to-severe inflammation are absent. These observations warrant further exploration in a broader IDA population.

Genome-wide Meta-analysis Identifies Novel Genes Associated with Recurrence and Progression in Non-muscle-invasive Bladder Cancer.

BACKGROUND: Non-muscle-invasive bladder cancer (NMIBC) is characterized by frequent recurrences and a risk of progression in stage and grade. Increased knowledge of underlying biological mechanisms is needed. OBJECTIVE: To identify single nucleotide polymorphisms (SNPs) associated with recurrence-free (RFS) and progression-free (PFS) survival in NMIBC. DESIGN, SETTING, AND PARTICIPANTS: We analyzed outcome data from 3400 newly diagnosed NMIBC patients from the Netherlands, the UK, Canada, and Spain. We generated genome-wide germline SNP data using Illumina OmniExpress and Infinium Global Screening Array in combination with genotype imputation. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: Cohort-specific genome-wide association studies (GWASs) for RFS and PFS were performed using a Cox proportional hazard model. Results were combined in a fixed-effect inverse-variance weighted meta-analysis. Candidate genes for the identified SNP associations were prioritized using functional annotation, gene-based analysis, expression quantitative trait locus analysis, and transcription factor binding site databases. Tumor expression levels of prioritized genes were tested for association with RFS and PFS in an independent NMIBC cohort. RESULTS AND LIMITATIONS: This meta-analysis revealed a genome-wide significant locus for RFS on chromosome 14 (lead SNP rs12885353, hazard ratio [HR] C vs T allele 1.55, 95% confidence interval [CI] 1.33-1.82, p = 4.0 × 10-8), containing genes G2E3 and SCFD1. Higher expression of SCFD1 was associated with increased RFS (HR 0.70, 95% CI 0.59-0.84, pFDR = 0.003). Twelve other loci were suggestively associated with RFS (p < 10-5), pointing toward 18 additional candidate genes. For PFS, ten loci showed suggestive evidence of association, indicating 36 candidate genes. Expression levels of ten of these genes were statistically significantly associated with PFS, of which four (IFT140, UBE2I, FAHD1, and NME3) showed directional consistency with our meta-analysis results and published literature. CONCLUSIONS: In this first prognostic GWAS in NMIBC, we identified several novel candidate loci and five genes that showed convincing associations with recurrence or progression. PATIENT SUMMARY: In this study, we searched for inherited DNA changes that affect the outcome of non-muscle-invasive bladder cancer (NMIBC). We identified several genes that are associated with disease recurrence and progression. The roles and mechanisms of these genes in NMIBC prognosis should be investigated in future studies.

Semi-Quantitative Multiplex Profiling of the Complement System Identifies Associations of Complement Proteins with Genetic Variants and Metabolites in Age-Related Macular Degeneration.

Age-related macular degeneration (AMD) is a major cause of vision loss among the elderly in the Western world. The complement system has been identified as one of the main AMD disease pathways. We performed a comprehensive expression analysis of 32 complement proteins in plasma samples of 255 AMD patients and 221 control individuals using mass spectrometry-based semi-quantitative multiplex profiling. We detected significant associations of complement protein levels with age, sex and body-mass index (BMI), and potential associations of C-reactive protein, factor H related-2 (FHR-2) and collectin-11 with AMD. In addition, we confirmed previously described associations and identified new associations of AMD variants with complement levels. New associations include increased C4 levels for rs181705462 at the C2/CFB locus, decreased vitronectin (VTN) levels for rs11080055 at the TMEM97/VTN locus and decreased factor I levels for rs10033900 at the CFI locus. Finally, we detected significant associations between AMD-associated metabolites and complement proteins in plasma. The most significant complement-metabolite associations included increased high density lipoprotein (HDL) subparticle levels with decreased C3, factor H (FH) and VTN levels. The results of our study indicate that demographic factors, genetic variants and circulating metabolites are associated with complement protein components. We suggest that these factors should be considered to design personalized treatment approaches and to increase the success of clinical trials targeting the complement system.

Integrating Metabolomics, Genomics, and Disease Pathways in Age-Related Macular Degeneration: The EYE-RISK Consortium.

PURPOSE: The current study aimed to identify metabolites associated with age-related macular degeneration (AMD) by performing the largest metabolome association analysis in AMD to date, as well as aiming to determine the effect of AMD-associated genetic variants on metabolite levels and investigate associations between the identified metabolites and activity of the complement system, one of the main AMD-associated disease pathways. DESIGN: Case-control association analysis of metabolomics data. PARTICIPANTS: Five European cohorts consisting of 2267 AMD patients and 4266 control participants. METHODS: Metabolomics was performed using a high-throughput proton nuclear magnetic resonance metabolomics platform, which allows quantification of 146 metabolite measurements and 79 derivative values. Metabolome-AMD associations were studied using univariate logistic regression analyses. The effect of 52 AMD-associated genetic variants on the identified metabolites was investigated using linear regression. In addition, associations between the identified metabolites and activity of the complement pathway (defined by the C3d-to-C3 ratio) were investigated using linear regression. MAIN OUTCOME MEASURES: Metabolites associated with AMD. RESULTS: We identified 60 metabolites that were associated significantly with AMD, including increased levels of large and extra-large high-density lipoprotein (HDL) subclasses and decreased levels of very low-density lipoprotein (VLDL), amino acids, and citrate. Of 52 AMD-associated genetic variants, 7 variants were associated significantly with 34 of the identified metabolites. The strongest associations were identified for genetic variants located in or near genes involved in lipid metabolism (ABCA1, CETP, APOE, and LIPC) with metabolites belonging to the large and extra-large HDL subclasses. Also, 57 of 60 metabolites were associated significantly with complement activation levels, independent of AMD status. Increased large and extra-large HDL levels and decreased VLDL and amino acid levels were associated with increased complement activation. CONCLUSIONS: Lipoprotein levels were associated with AMD-associated genetic variants, whereas decreased essential amino acids may point to nutritional deficiencies in AMD. We observed strong associations between the vast majority of the AMD-associated metabolites and systemic complement activation levels, independent of AMD status. This may indicate biological interactions between the main AMD disease pathways and suggests that multiple pathways may need to be targeted simultaneously for successful treatment of AMD.

Sex-Specific Regulation of Inflammation and Metabolic Syndrome in Obesity.

OBJECTIVE: Metabolic dysregulation and inflammation are important consequences of obesity and impact susceptibility to cardiovascular disease. Anti-inflammatory therapy in cardiovascular disease is being developed under the assumption that inflammatory pathways are identical in women and men, but it is not known if this is indeed the case. In this study, we assessed the sex-specific relation between inflammation and metabolic dysregulation in obesity. Approach and Results: Three hundred two individuals were included, half with a BMI 27 to 30 kg/m2 and half with a BMI>30 kg/m2, 45% were women. The presence of metabolic syndrome was assessed according to the National Cholesterol Education Program-ATPIII criteria, and inflammation was studied using circulating markers of inflammation, cell counts, and ex vivo cytokine production capacity of isolated immune cells. Additionally, lipidomic and metabolomic data were gathered, and subcutaneous fat biopsies were histologically assessed. Metabolic syndrome is associated with an increased inflammatory profile that profoundly differs between women and men: women with metabolic syndrome show a lower concentration of the anti-inflammatory adiponectin, whereas men show increased levels of several pro-inflammatory markers such as IL (interleukin)-6 and leptin. Adipose tissue inflammation showed similar sex-specific associations with these markers. Peripheral blood mononuclear cells isolated from men, but not women, with metabolic syndrome display enhanced cytokine production capacity. CONCLUSIONS: We identified sex-specific pathways that influence inflammation in obesity. Excessive production of proinflammatory cytokines was observed in men with metabolic syndrome. In contrast, women typically showed reduced levels of the anti-inflammatory adipokine adiponectin. These different mechanisms of inflammatory dysregulation between women and men with obesity argue for sex-specific therapeutic strategies.

Measurement and genetic architecture of lifetime depression in the Netherlands as assessed by LIDAS (Lifetime Depression Assessment Self-report).

BACKGROUND: Major depressive disorder (MDD) is a common mood disorder, with a heritability of around 34%. Molecular genetic studies made significant progress and identified genetic markers associated with the risk of MDD; however, progress is slowed down by substantial heterogeneity as MDD is assessed differently across international cohorts. Here, we used a standardized online approach to measure MDD in multiple cohorts in the Netherlands and evaluated whether this approach can be used in epidemiological and genetic association studies of depression. METHODS: Within the Biobank Netherlands Internet Collaboration (BIONIC) project, we collected MDD data in eight cohorts involving 31 936 participants, using the online Lifetime Depression Assessment Self-report (LIDAS), and estimated the prevalence of current and lifetime MDD in 22 623 unrelated individuals. In a large Netherlands Twin Register (NTR) twin-family dataset (n ≈ 18 000), we estimated the heritability of MDD, and the prediction of MDD in a subset (n = 4782) through Polygenic Risk Score (PRS). RESULTS: Estimates of current and lifetime MDD prevalence were 6.7% and 18.1%, respectively, in line with population estimates based on validated psychiatric interviews. In the NTR heritability estimates were 0.34/0.30 (s.e. = 0.02/0.02) for current/lifetime MDD, respectively, showing that the LIDAS gives similar heritability rates for MDD as reported in the literature. The PRS predicted risk of MDD (OR 1.23, 95% CI 1.15-1.32, R2 = 1.47%). CONCLUSIONS: By assessing MDD status in the Netherlands using the LIDAS instrument, we were able to confirm previously reported MDD prevalence and heritability estimates, which suggests that this instrument can be used in epidemiological and genetic association studies of depression.

Identification of ADHD risk genes in extended pedigrees by combining linkage analysis and whole-exome sequencing.

Attention-deficit/hyperactivity disorder (ADHD) is a common neurodevelopmental disorder with a complex genetic background, hampering identification of underlying genetic risk factors. We hypothesized that combining linkage analysis and whole-exome sequencing (WES) in multi-generation pedigrees with multiple affected individuals can point toward novel ADHD genes. Three families with multiple ADHD-affected members (Ntotal = 70) and apparent dominant inheritance pattern were included in this study. Genotyping was performed in 37 family members, and WES was additionally carried out in 10 of those. Linkage analysis was performed using multi-point analysis in Superlink Online SNP 1.1. From prioritized linkage regions with a LOD score ≥ 2, a total of 24 genes harboring rare variants were selected. Those genes were taken forward and were jointly analyzed in gene-set analyses of exome-chip data using the MAGMA software in an independent sample of patients with persistent ADHD and healthy controls (N = 9365). The gene-set including all 24 genes together, and particularly the gene-set from one of the three families (12 genes), were significantly associated with persistent ADHD in this sample. Among the latter, gene-wide analysis for the AAED1 gene reached significance. A rare variant (rs151326868) within AAED1 segregated with ADHD in one of the families. The analytic strategy followed here is an effective approach for identifying novel ADHD risk genes. Additionally, this study suggests that both rare and more frequent variants in multiple genes act together in contributing to ADHD risk, even in individual multi-case families.